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Abstract 

We apply our recently developed resonance perturbation theory to describe the 
dynamics of magnetization in paramagnetic spin systems interacting simultane- 
ously with local and collective bosonic environments. We derive explicit expressions 
for the evolution of the reduced density matrix elements. This allows us to calculate 
explicitly the dynamics of the macroscopic magnetization, including characteristic 
relaxation and dephasing time-scales. We demonstrate that collective effects (i) 
do not influence the character of the relaxation processes but merely renormalize 
the relaxation times, and (ii) significantly modify the dephasing times, leading 
in some cases to a complicated (time inhomogeneous) dynamics of the transverse 
magnetization, governed by an effective time-dependent magnetic field. 

1 Introduction 

When quantum systems interact with their environments the effects of relaxation and 
decoherence occur [1-8]. In this paper we study relaxation and decoherence in quan- 
tum macroscopic systems of "effective" spins interacting simultaneously with both local 
and collective thermal environments. By "effective" we mean that our approach can 
be applied not only to magnetic spin systems, but also to many quantum systems with 
discrete energy levels, including recently widely discussed quantum bits (qubits) based 
on superconducting Josephson junctions and SQUIDs [9-14]. We also would like to 
mention here the research on ephaptic coupling of cortical neurons, when both local 
and collective electrical fields play a significant role in the synchronization dynamics 
of neurons [15]. We assume that spins do not interact directly among themselves, but 
only through their interactions with collective (energy conserving and energy exchange) 
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bosonic environments ( "thermal baths" , "reservoirs" ) . The relaxation in these systems 
is caused by energy exchange between the environments and spins. The rate of relax- 
ation is usually characterized by the spectral density of noise of the reservoirs at the 
transition frequency, uj, of spins in their effective magnetic field, and by the interac- 
tion constant between a spin and an environment [7,9,14]. The rate of decoherence 
usually has a more complicated dependence on the parameters of spins, their local en- 
vironments, and interaction constants [7-9,14,16]. In particular, low-frequency noise 
(1// noise) makes a significant contribution to the decoherence rate [17-20] (see also 
references therein). 

Usually relaxation and decoherence are unwanted effects, for example, in a quantum 
computer one must maintain quantum coherence for long times [5]. But dissipative 
effects can also be put to good use, for example, in magnetic resonance imaging (MRI) 
[21-23]. Indeed, in this case, different values of relaxation times (Ti) for different 
substances (e.g. water and biological tissues) allow one to distinguish and visualize 
pathological developments in tissues [21]. Dissipative effects can also be utilized, for 
example, to analyze and classify the influence of many types of defects and impurities, 
in order to improve the properties of materials [24] . 

Improving our understanding of relaxation and decoherence processes is important 
for many fields of science and for many applications. The main problem associated 
with dissipative effects is that there are many different sources of noise and thermal 
fluctuations which lead to relaxation and decoherence. We mention only some of them, 
electromagnetic and acoustic fliictuations (bosonic degrees of freedom), magnetic fluc- 
tuations (such as two- level systems in superconducting materials), charge defects, and 
non-equilibrium quasiparticles. Generally, it is impossible to eliminate all sources of 
noise, so some additional classification can be useful. For example, in [16] we demon- 
strated that, for a system of N spins interacting with bosonic environments, one can 
introduce clusters of reduced density matrix elements in such a way that to a given 
cluster corresponds a decoherence rate describing the fading of all matrix elements 
belonging to it. When dealing with a quantum algorithm in quantum computation, 
this could imply that the decay of some clusters is rapid, but - if the algorithm is 
built mainly on the use of slower decaying clusters - that decay may not influence 
significantly the fidelity of the quantum protocol. 

In this paper we are mainly interested in the effects produced by simultaneous 
influence of both local and collective bosonic environments on the dynamics of a col- 
lective magnetization in a system of N non-interacting (paramagnetic) spins in a time- 
independent magnetic field. The local and collective environments include both energy 
conserving and energy exchange interactions with spins. This allows us to determine 
conditions of applicability of the Bloch equation for describing the evolution of the 
magnetization. We also consider two (and more) ensembles of spins with different pa- 
rameters and strengths of interactions with their environments. Using our approach 
based on resonant perturbation theory [16], we derive explicit expressions for the time 
evolution of the reduced density matrix elements and, consequently, for the macro- 
scopic magnetization. We explicitly calculate the relevant relaxation and decoherence 
rates. The obtained results are important for many applications including MRI and 
for studying collective effects in materials for superconducting qubits. 
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Main results of the paper 



• Single spin dynamics. We consider a microscopic, Hamiltonian model of N spins 
interacting with local and collective bosonic thermal reservoirs, via energy conserving 
and energy exchange interactions. In Theorem 2.1 wc derive a rigorous expression 
for the reduced density matrix of a single spin, consisting of a main term describ- 
ing relaxation and dephasing, plus a remainder term which is small in the couplings 
homogeneously in time. 

• Single spin relaxation. We show that the single-spin relaxation rate is given by 

7reiax = ^ coth(/3a;/2) {X^Jg^{oj) + ii^Jg^{uj)] , 

where uo is the spin frequency, A and /v, arc the strengths of the energy exchange collective 
and local couplings, respectively, and where Jg{u) is the reservoir spectral density. Only 
energy-exchange couplings contribute to this rate, and the effect of the local and the 
collective reservoirs are the same. 

• Single spin dephasing. We show that the single-spin dephasing rate is given by 

-1 

Tdeph — ^"^Telscyi + Tcons + T ) 

where 7cons is a contribution stemming only from the energy conserving local and collec- 
tive interactions, determined by the spectral density of the reservoir at zero frequency 
(see (3.3)). The contribution 7' encodes the effect on dephasing of a single spin due 
to all other spins. It is defined as follows. The time-dependence of the single spin 
off-diagonal density matrix elements has a very complicated, not exponentially decay- 
ing contribution coming from the collective coupling. The term 7' is defined to be the 
reciprocal of the time by which that quantity is reduced to half its initial value. 

The explicit expression of 7' is not simple (see (2.33), (2.31)). For small ratio 
r between the strengths of the collective to the local couplings we have 7' = O(r^) 
(independent of the number N of spins). For large collective coupling we have 7' ~ 
const. 7reiax) for a constant not depending on N . 

• Evolution of magnetization. We consider the spins in a homogeneous magnetic 
field pointing in the z-direction. We show that the z-component of the total magne- 
tization vector relaxes to its equilibrium value at the single -spin relaxation rate 7rclax* 
This verifies the correctness of the usual Bloch equation (3.13) for the z-component. 
The equation for the transverse total magnetic field is given by a modified Bloch equa- 
tion (3.15), with a time-dependent dephasing time (T2 = 72 (t)) and a time-dependent 
effective magnetic field. For large times, the coefficients in the modified Bloch equation 
approach stationary values and give rise to the usual Bloch equation with renormalized 
r2(oo) time and renormalized effective magnetic field. We show that 

7777 \" ~ oTrelax + Tcons + (-^ ~ 1)7 ) 
i2(OOj 2 

where 7" > is independent of N . For small ratio r between the strengths of the col- 
lective to the local couplings we have 7" = O(r^). Consequently, if r ~ N~^/^ then the 
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collective coupling gives a non-vanishing renormalization to the (asymptotic) T2 time, 
while if r ~ jV"^/^^'^ (any e > 0) or smaller, then no collective effect is visible in the de- 
phasing. An interesting question is what happens for r ~ N~^/^'^^ or larger. Then the 
expression for T2(oo) suggests that the collective interaction may decrease the T2 time 
drastically for large N. However, this range of interaction parameters is not accessible 
by our perturbation theory approach, and more work in this direction is required. It 
is important to note here that in order to derive our rigorous result, Theorem 2.1, we 
need a strong smallness condition on all coupling constants (see (2.11)). As explained 
in Section 5, wc expect that our result should hold for collective coupling constants 
up to size 0(A'"~^/^) and local coupling constants of size 0{N^) (relative to the spin 
frequency). However, for r = 0{N~^/^'^^) we do not think that usual perturbation 
theory can be applied, and a different approach should be taken. 

Wc also examine the situation where wc have two (or more) species of spins, A and 
B, each species coupled homogeneously to local and collective reservoirs (with a single 
collective reservoir for both species). We show that the z-component of the magne- 
tization of either species relaxes with single-spin relaxation time (associated to that 
species). The transverse magnetization dephases following a modified Bloch equation 
with time-dependent T2-time and effective magnetic field. For large times, the r2-time 
of species A approaches the limiting value 

7f-^7 r = ;^7relax,A + 7cons,A + {Na - 1)ja + NbJB, 

J^2,a[oo) z 

where Na and Nb are the number of spins in each class, and jai 7b > 0. For small ratio 
rA, fB of the collective and local coupling constants, we have 7 a = O(r^), 75 = 0{r^). 
The total magnetization is the sum of that of species A and B. It is the sum of two 
terms decaying (relaxing and dephasing) at different rates, and so we cannot associate 
to it a total relaxation time or a total dephasing time. 

The effects of collective interactions between effective spins and thermal environ- 
ments, discussed in this paper, can represent a significant interest, for example, in 
NMR, MRI, and quantum computation. The presence of energy conserving and energy 
exchange collective effects can be investigated experimentally, for example, in NMR 
experiments by (i) creation and controlling of collective effects and (ii) analyzing relax- 
ation and dephasing time-scales and time-dependencies of magnetization as functions 
of characteristic parameters. 
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2 Model, single spin dynamics 

We consider N non interacting spins 1/2 coupled to local and collective bosonic heat 
reservoirs. The full Hamiltonian is given by 

N N 

H = -h^UnS^ + Y^HKr^+Hn (2.1) 

n=l n=l 

N N 

+ J2>^nS^^M9c) + J2''-^n^Mfc) (2.2) 
n=l n=l 
N N 

+ J2 l^nSl ® Man) + ""^^n ® Mfn)- (2-3) 

n=l n=l 

Below wc use dimensionless variables and parameters. To do so, we introduce a 
characteristic frequency, ljq, typically of the order of spin transition frequency. The 
total Hamiltonian, energies of spin states, and temperature are measured in units Huq. 
The frequencies of spins, a;„ > 0, bosonic excitations, oj{k) = c\k\ (where c is the 
speed of light) , the wave vectors of bosinic excitations are normalized by loq /c, and all 
constants of interactions are measured in units coq. A dimemsionless time is defined as 

t — )• UJot. 

In (2.2), (2.3), a;„ > is the frequency of spin n, 



1 
-1 



and S'' = ]- 
2 



1 

1 



(2.4) 



and Sn^ denotes the S^'^ of spin n. is the Hamiltonian of the bosonic collective 
reservoir, 

Hn= [ \k\a*ik)aik)d^k, (2.5) 

and i^R„ is that same Hamiltonian pertaining to the n-th individual reservoir. For a 
square-integrable form factor h{k), A; G M^, 0(/i) is given by 

(f){h) = ^ [ {h{k)a*{k) + h{k)*a{k)}d^k. (2.6) 
V2 Jm.3 

The real numbers A„, x„, f„ are coupling constants, measuring the strengths of 
the various interactions as follows: 

Xn energy exchange collective coupling 

x„ energy conserving collective coupling 

Hn energy exchange local coupling 

Vn energy conserving local coupling 

We introduce the maximal size of all couplings, 

a := max{|x„|, |A„|, |z/n|}. (2.7) 
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The energies of the N uncoupled spins are the eigenvalues of = — J2n=i ^nSn, 
given by -i J2n=i ^n^n, where o-„ G {1, -1}. We denote hy Lpa_ = i^a^® ■ ■ ■ ® ^aj^ the 
corresponding eigenvector. Bohr energies (energy differences) are thus given by 



1 ^ 

e(£,r) = --^w„(c7„-T„). (2.8) 

n=l 

Assumptions. 

(A) We consider the spin frequencies {oJn} to be uncorrclatcd in the following sense: 

-(f e(a,T) = e{^,¥) then cr„ — r„ = — t!^ for all n. (2-9) 

In particular, we do not allow any of the cj„ to be the equal. However, we can 
describe a homogeneous magnetic field within the constraint (2.9) by considering 
a distribution a;„ = w + for some fluctuation having, say, uniform distri- 
bution in some interval. Then assumption (A) is satisfied almost surely. In view 
of such fluctuations, relation (2.9) is reasonable from a physical point of view, and 
its mathematical advantage is that it breaks permutation symmetry and hence 
reduces the degeneracies of the energies e. 

(B) The smallest gap between different Bohr energies of the non-interacting spin 
energies (2.8) is 



^ 1 . 

A = — mm 
2 mn,m'„ 



N 



\{0}, (2.10) 



where the minimum is taken over sequences G {~2, 0, 2}. As our approach 

is based on perturbation theory of Bohr energy differences, their displacement 
under interaction, which is of size N'^a^, should be small relative to A, 

iV^A^ « A. (2.11) 

For a;„ = w constant, we have A = huj. Hence for a homogeneous magnetic field 
ujn = u + 6ujn with Small fluctuation {5iOn)/co « 1, we have A = lo + 0{{6uj)), and 
the r.h.s. of (2.11) is independent of N. Condition (2.11) is a serious restriction on 
the coupling strength for large systems (big N) . Our analysis uses this condition 
in several technical estimates of remainder terms, stemming from perturbation 
theory (see also [16]). However, it is seen from physical considerations, presented 
in Section 5, that the true condition should read a^N « u and a£ « co, where 
Oc and a£ are the sizes of collective and local coupling constants, and co is the 
typical frequency of a spin. 

(C) Regularity of form factors: denote by h any of the functions fc, dc, fn, dn in the 
Hamiltonian H. Let r > 0, S G S"^ be the spherical coordinates of M^. Then 
/i(r, S) = r^e~'^''"/i'(S), with p = — 1/2-1- n, n = 0,1,2,... and m = 1,2, and 
where h' is any angular function. (Less restrictive requirements on h are necessary 
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only [16], but they are more technical to describe, so we restrict our attention to 
h satisfying this condition. This family of form factors contains the usual physical 
ones [8].) 

Given e, (2.8), the number 

7Vo(e) = {n : (Tn = Tn for any {a,T) with e(a, r) = e} (2-12) 

depends on e alone, and the number of different configurations {a, r) with constant value 
(2.8) is 2^°^^\ There are 2^0(6) elements {'^a^Ps^x) of the (reduced) density matrix of 
the spins with fixed value (2.8). As shown in [16], these elements evolve in time jointly, 
and independently of elements associated with any other value of (2.8). 
We consider unentangled initial states 

P0 = PSi® ■ ■ ■ ® PSn® PKi® ■ ■ ■ ® PKn ® PR> 

where ps^ are arbitrary single spin states, and pR^. , /9r are thermal equilibrium states 
of single reservoirs, all at temperature T = 1/^ > 0. 
The reduced density matrix pf^ of spin j is given by 

the trace being taken over all spins j and over all reservoirs. 

Let Aj be an observable of the j-th spin, and denote its dynamics by 

(A,), = TVp(^)^,-, 

where the trace is taken over the space of Sj. Our goal is to find a representation of 

{Aj)^. For a square integrable form factor h{k) = h{\k\,T,) (spherical coordinates of 
M^), the spectral density of the reservoir associated to h is given by 

Jh{uj)=7ruj'^ [ |/i(w,E)|2dS. ^ (2.13) 

Decay rates are given by coupling constants squared times Jh{oj) coth(/3u;/2) at values 
uj corresponding to Bohr frequencies of the spin system. Energy-conserving processes 
are associated with the Bohr frequency a; = 0, and since coth(/3a;/2) ~ Ld~^ as w ~ 0, 
we introduce 

MO) = hm (2.14) 



^Let Chit) = i[((/>(/i)e"^a0(/i)e"'*^R)/j + {e'*^R<?!)(ft)e""^R(/>(/i));3] be the symmetrized correlation 
function of a reservoir in thermal equilibrium at temperature T = 1//3, with Hn and <f){h) given in (2.5) 
and (2.6). The Fourier transform Ch{w) = /q°° e~"^*C(t)dt, w > 0, is related to the spectral density by 

ReChioj) = Jh(w)coth(^w/2). 
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Define the quantities 

1 e/^^J 





- i^0:j—i{^p9M) + i^p9M)} (2-15) 

Cj = e-^'^^ (2.16) 

Z^j = Q-^^il-^ j^^fi^il-^ (2.17) 

-!-P.V. / ' ^ ' ' coth(;9|n|/2)dtx (2.18) 

= ^{AHc('^.) + /^''^P.(^i)}coth(;3a;,/2) + ^{x|j/,(0)+z.|j/^^ (2.19) 
With this notation in place wc have the foUowing result. 

Theorem 2.1 (Dynamics of single spin) For any observable Aj of spin j, j = 
1, . . . ,N , and t > 0, we have 

{Aj\ = Z^^.TTe-^'''^Aj (2.20) 

+e-%('^^+^) |[pi^-)]n - ^z^^} mn - [Ajh2) (2.21) 

^gi*(-.,+x,+iy,)^^.(^^^) [pi^'^]2i [A,]u (2.22) 

^^ii(_^^,+x,+iy,)c.(iv,i) [p<^^% [Aj]2i (2.23) 

+0(a2). (2.24) 



The quantity Cj{N,t) involves the interaction parameters and the initial condition of 
all spins other than j, 

CjiN,t) = n|[e^*-<+-e«^rll±^(a, + [pW]n(l-a^))+e"^r| (2.25) 
z± = i |i6/(l + ci) ± sj-b]il + q)2 + 4ai[ai - Ml - q)]| (2.26) 

ai = l + i-4 ^ 2.27 

In (2.26), the square root is the principal value (cut on the negative real axis), 6;, q are 
given in (2.15), (2.16), and 

a, = -ixfP.V. / ^M^dV (2.28) 
^ Jr3 \p\ 

Discussion of the factor Cj{N,t). 

Clearly Cj{N,0) = 1. For vanishing energy-conserving collective coupling, = (all 
I), we have Cj{N, t) = 1 for all t >0. This follows from a; = 0, z'^ = i6;(l + q), = 
and ai = — 1/q. 
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As soon as the collective energy-conserving coupling is switched on, ki ^ 0, the 

analysis of Cj{N,t) is difficult. The factors in the product (2.25) decay with rate at 
least 

^j=mm{Qz+,Qz-}, (2.29) 

so for a homogeneous system (each factor the same) we have the estimate \C{N,t)\ < 
Cive-T(^-i)*, with 

7 = min7j. (2.30) 

j 

Of course, Cj does not depend on j anymore. This estimate says that C decays in time 
with rate 7(A'' — 1), but we have a prefactor depending on N. We have the upper bound 

Cn < e(^-i)<=', for 



c' = In <^ 2 



'^"■(a + [yOo]ii(l - a] 



+ 1 ^ > 0. (2.31) 



1 + ca^ 

We know that for t = the true upper bound on \C{N, t)\ corresponds to Cjv = 1, but 
this does not mean at all that \C{N,t)\ < e~'''(^~^)*. The estimate 

|C(iV,i)| < e(^-i)[-T*+'^'l, with 7 = min^ ^^-, see (2.29) (2.32) 

shows that \C\ decays to 1/2 (half of its initial value) no later than at time J^^[^^+c']. 
We thus call 



7' = 7 



In 2 ,1"^ 
+ c' 



^ _ ^ _ (7' «^ 7/c' for large N) (2.33) 

the decay rate of of \C{N,t)\. Let us examine this decay rate in the two cases where 

is either very small or very close to one. (2.34) very small corresponds to the situation 
where the collective interactions (x ^ A) are much smaller than the local ones {11 ^ v). 
The situation where (2.34) is unity describes very large collective coupling relative to 
the local ones. 

For small collective coupling, r « 0, we obtain 7' ~ const. rx, where const, does not 
depend on N. In the limit r ^ we get 7' = 0, which is the correct behaviour as we 
have seen above (C = 1 in this setting, no decay). 

For large collective coupling, r « 1, we obtain 7' ^ const. 6, where const, does not 
depend on A*" and b is given in (2.15). 



Comparison to exactly solvable model. 

If the spins interact with the reservoirs only through energy conserving channels, then 
-^n = A^n = ill (2-2), (2.3). This model is exactly solvable. By proceeding as in [16] 
( "Resonance theory of decoherence and thermalization" , proof of Proposition 7.4) , one 
finds the following exact formula for the evolution of the reduced density matrix element 
of a single spin. For simplicity of notation, we take all x„ to be constant Xc (collective) 
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and all i/„ constant ui (local). We also take all local form factors equal (J^) and all 
collective ones too (/c). Then we find 

N 

The sum is over ak = ±1/2, k = l,...,N,k^j, where aj = 1/2 corresponds to the 
energy eigenstate ipi = [1 0]-^ of (see (2.4)). The decoherence functions and Lamb 
shift are given by 

m = [ |/(fc)pcoth(/3|fc|/2) ^^"'|iy^ d^fc (2.36) 

Jr3 |fe|^ 

S(t) = (2.37) 

Of course, the populations are time-independent in this model, [pp^];^ = [Po'^jii for all 
t > 0, I = 1,2. The time-dependence in the exponentials in (2.35) becomes linear for 
large times, T{t) — ^ tJ{0) and >c'^S{t) — ta as t — )■ oo, where J(0) and a are given in 
(2.14) and (2.28) with replaced by Xc (see [16]). 

It is not hard to see that upon the replacements T{t) i-)- tJ{0) and >c1S{t) i— )■ ta the 
exact formula (2.35) coincides precisely with expression (2.20)-(2.23) for A s.t. = 
[A]22 = [A]2i = 0, [A]i2 = 1 (so that {A)t = [pthi)- The factor Cj(iV,t) = C{N,t) is 
thus identified with the sum of the product in (2.35). If all spins are initially in the 
same state, characterized by the population probability < p < 1 for the state with 
a = 1/2, we obtain 

C{N, t) = [pe-'"* + (1 - p^"^] , with a given in (2.28). 

Clearly |C(A^, t)| < 1 for all times and all A'^. Also, for all n G Z, we have \C{N, n7ra~^)| = 
1 and |C(iV, (n + ^)7ra-^)| = |1 - 2p|^-i for iV large and p / 0, 1. Therefore the 
factor C{N,t) oscillates in size between zero and one, with frequency \a\/Tr proportional 
to the square of the energy-conserving collective coupling x^. 

3 Evolution of single spins and of magnetization 
3.1 Single spin releixation and dephasing times 

The term on the r.h.s. of (2.20) is the equilibrium average at temperature T = 
From (2.21) we obtain the relaxation rate of spin j, namely 7reiaxj = bj{cj + 1). The 
single spin relaxation rate is 

Trelaxj = l/7"rclax,i = ^ COth(/3Wj/2) {\^jjg^{uj) + /i^ J^^ {uj)] . (3.1) 

The single-spin relaxation time depends on the local {^ij) and collective (Aj) couplings 
in the same manner: In the relaxation process, the collective reservoir acts as a local 
reservoir. 
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Next we consider the dcphasing time determined by (2.22), (2.23). There are two 
contributions to the time decay. One comes from spin j itself and is given by Yj , the 
other one comes from all other spins than j and is given by Cj. One sees from (2.26) 
that Qzf" > 0, and that minj^Jz;"*", ^Jz^"} = aibi = 0. It follows that if the energy 
conserving collective coupling and at least one of the energy-exchange couplings (local 
or collective) do not vanish (so that aibi / 0), then we have (2.32) with 7 > 0. The 
single-spin dephasing rate is thus Yj + 7' which we can write as 

TdephJ — ^IrelaxJ ~^ Tconsj "I" T ) C^-^) 

where 

Tconsj = ^{^jJfM + '^jJfM (3.3) 

is a contribution coming purely from the energy-conserving interactions, in which the 
local and collective couplings play the same role. The last term in expression (3.2) is 
due to the presence of the TV — 1 spins other than the considered one. As we have 
seen after (2.34), if the collective coupling is small (r ~ 0), then 7' ~ x^r « and 
hence the last term in (3.2) is negligible. If the collective coupling is large (r ~ 1), then 

y ~ 6 ~ TrelaxJ- 

Conclusions. • The single-spin relaxation rate is the sum of two contributions 
from the local and the collective energy-exchange interactions (3.1). The collective 
term has the same form as the local term, and the presence of all other spins does not 
influence the single spin relaxation rate. 

• The single-spin dcphasing rate has three contributions (3.2). One is half the 
relaxation rate (exchange interactions), one comes from energy conserving interactions 
(local and collective), and a third term which is due to the presence of all other spins. 
That last term (7') is negligible for small collective coupling, and renormalizes the 
dephasing rate for strong collective couplings by an amount independent of the number 
of spins. 

3.2 Evolution of magnetization 

Let 



S = 



Sy 



be the total magnetization vector, where 5^'^'^ = '^jLi Sj'^'^. It is convenient to 
introduce the complex (non-hermitian) observable 

We use Theorem 2.1 with Aj = S^'^'"" to obtain 

= ^ tanh(;3a;^V2)[l - e-*/^-'--^] + e"*/^-'-'^' (S^)^ + 0{a^), (3.4) 
{S-)t = e'*(-'^^+^^+'^^)C,(Ar,t)(5-)o + 0(a2)_ (35) 
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Purely local coupling. In the absence of collective coupling (A„ = = x„), the 
above equations simplify to 

= ^tanh(;9a;j-/2)[l-e-*/^-'--^] + e-*/^-'-.^(5|>Q + 0(a2), (3.6) 
{Sr)t = e'*(-'^^+^^+'^^)(5-)o + 0(a2). (3.7) 

where Treiaxj is given by l/7rciaxjr (3.1) with Xj = 0, Xj,Yj are given in (2.18), (2.19) 
with Xj = jjLj = 0. The factor Cj{N,t) equals 1 (as discussed after (2.28)). 

3.2.1 Homogeneous magnetic field 

In this section we derive the evolution of the magentization vector in a homogeneous 
magnetic field, characterized by luj = to + Sojj with Suj — >■ (see also assumption (B) 
after (2.11)). This is the description of an elementary volume of many spins sitting in 
a magnetic field with gradient much smaller than the size of the elementary volume. 

We consider all spins initially in the same state. We take all local couplings to be 
the same, i.e., gj = gg and all collective couplings to be the same, fj = fc, and all 
coupling constants independent of j. In this limit, we have in formulas (3.4), (3.5) 

ojj = CO, Xj = X, Yj = 

where Trdax, X, Y arc given in (3.1), (2.18), (2.19), with ujj,fj,gj and all coupling 
constants replaced by their constant values, in particular, 

7rels« = = \ COth(^a;/2) {X^ JgS^^) + /X^Jg,(w)} • (3.8) 

Furthermore, we have Cj{N,t) = C{N,t), with (see (2.25)) 

C{N,t) = [V{t)]^-^ (3.9) 
V{t) = U^+-e«^"l :^^^(Q + [po]ii(l-a))+e'*^". (3.10) 

We sum equations (3.4) and (3.5) over j to obtain (dropping the O(a^) terms) 

{S')t = y tanh(/3L^/2)[l -e-*/^-i-] +e-*/""-i- (5^)o (3.11) 
{S-)t = e'*(-'"+^+'^)[P(t)]^-i(5-)o. (3.12) 
It is clear that (3.11) is the integrated version of the Bloch equation 

A (S^)^ = [ {S]), - y tanh(^a;/2)] (3.13) 

corresponding to the homogeneous magnetic field B = B^Cz = —^&zi with relaxation 
time 

^1 ~ "^relax ~ 1 / Trelax ) 
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see (3.8). The Bloch equation for the transverse magnetization would read 

^(5-), = -^(5-), + ii?,(5-),. (3.14) 

However the true evolution, (3.12), is not of this form. By differentiating (3.12) we 
obtain 

= -m{S-)t + iB{t){S-)t, (3.15) 

with 

r(t) = ^7rele« + 7cons-(iV-l)Re^lnP(t), (3.16) 

B{t) = -uj + X + {N -l)lm^lnV{t), (3.17) 

where 7rciax is the single-spin relaxation rate (3.1) and 7cons is the single-spin dephasing 
rate due to the energy-conserving interactions (3.3). 

Comparing (3.15) with (3.14) leads us to the identification of a time- dependent 
dephasing time T2 = l/r(t) and a time- dependent effective magnetic field Bz = B{t). 

The deviation of the true equation of evolution from the Bloch equation is given by 
the terms d/dtlnP(t) in (3.16), (3.17). We now estimate the size of this term for weak 
collective coupling, where r is small, see (2.34). It is not hard to see that 



< C\r\, lim 5^ = iz- = 4i6r tanh(^a;/2) + 0{r^), (3.18) 
t—^co 'D{t) 



m 

for a constant C independent of t (and N). 

Conclusions. The Bloch equation for the total magnetization (homogeneous 
magnetic field) holds with relaxation time Ti given by the single-spin relaxation rate 
(3.1) (no influence of the other spins). The total magnetization dephases with a time- 
dependent r2-time, T2 = T2{t). We have 1/T2(t) = i7reiax + 7cons + (A^- l)Rey'(t), see 
also (3.1), (3.3). The time-dependent part Y'{t) stems from the collective interaction. 
For weak collective interaction, r small (see (2.34)), we have |l^'(t)| < C\r\ (all times). 

The term {N - l)Re^P(t) in (3.16) is 0{Nr'^) for large times, see (3.18). If r is 
of the order l/y/N then this is of order one, and the collective interaction gives time- 
dependent modification of the dephasing time T2 with an asymptotically renormalized 
value ^ 

l/T2(oo) = -7rei£« + 7cons + {N - l)lmz- , Imz' = 0{r^). 

If r is smaller than N~^~'' (any e > 0) then the collective interaction has no effect 
in (3.16), (3.17) and the Bloch equation for transversal relaxation holds with T2 the 
single-spin dephasing time [7reiax/2 + 7cons]~^- For larger collective interaction we may 
get large corrections to the Bloch equation, since the last terms in (3.16), (3.17) may 
become large (big N). This regime does not enter the present perturbative setup, and 
more work on this issue is needed. 

Note that in any event, since lmz~ > 0, the collective interactions can only accel- 
erate the dephasing process. 
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3.2.2 Multi-species inhomogeneity 

Consider the situation where A'^ spins are grouped into two (or more) classes A and 
B. We describe the situation where within each class, the spins are homogeneous. We 
have two magnetic fields loa, wb, two sets of coupling constants (A^, As etc), two sets 
of form factors {gc,A, gc,B, QA, 9b etc). Let Na and Nb be the relative sizes, 

Na + Nb = N. 

If spin j belongs to class ^4, then (2.25) becomes 

Cj{N,t) = [VA{t)f^-' [VB{t)f'' , (3.19) 

with T>A{t),T>B{t) given as in (3.10) for species A^B. Let 

Sa= 5Z 

j in class A 

and, correspondingly, for the three components of this vector. We sum (3.4) and (3.5) 
over all indices of spins belonging to class A to obtain 

{S\)t = ^ tanh(/3wA/2)[l - c"*/^-'-.^] + e-*/^-'-.^(51)o + 0{a^) (3.20) 

(-S^)t = ^'^-''^+''^+''^^\VA{t)f^-\VB{t)f^Sj)^ + 0{a^). (3.21) 

Hence class A relaxes with single-spin relaxation time T^e\ax.,A according to the usual 
Bloch equation (3.18). For the transverse magnetization we obtain again a modified 
Bloch equation with time-dependent relaxation time and effective magnetic field, 

= -TAmSA)t + ^BA{t){S^)t + 0{a% (3.22) 

with 

^Ait) = ^7reiax,A + 7cons,A-(A^A-l)Re^lnpA(i)-A^BRe^lnPs(t) (3.23) 

BA{t) = -0JA + XA + {NA-l)lm^^lnVA{t)+NBlin^VB{t). (3.24) 

As in the previous paragraph, we see that for weak collective coupling and large times, 
TAit) converges to ^7reiax,A + 7cons,A - {Na - l)Imz^, and BA(t) converges to -uja + 
Xa + {Na — l)Imz^. We thus obtain the (asymptotic) dephasing rate for species A, 

7deph,A(oo) = ^7relax,^ + 7cons,^ + {^A " l)l^Z^ + NbI^z]^. ^ (3.25) 

straightforward generalization to s species Ai, . . . , As with sizes Nai -|- • • • -|- Na^ = N gives the 
transverse relaxation (i.e., dephasing) rates 

7deph,A. (OO) = + {Naj - l)lmzAj + ^NA^lrazAk, 

for j = 1, . . . , s. We conclude that the relaxation rate of each species is a single-spin relaxation rate, 
while the dephasing contains collective effects. In particular, the dephasing rate of class Aj depends 
on all other classes. 



14 



Recall again that for small collective interaction, ImZj^^ = 0{r\^), (2.34). 

Conclusions. The z-component of the total magnetization of each species A 
and B evolves according to the Bloch equation (3.13) with single-spin relaxation rates 

7relax,A and 7relax,B, (3.1). 

The transverse total magnetization of species A evolves according to a modified 

Bloch equation (3.22) (similarly for B). The dcphasing time becomes time-dependent 
(3.23), and takes the value T2^^(oo) = l/7deph,A(oo), (3.25) for large times and small 
collective coupling. 

The total magnetization is the sum of that of species A and B, {S)t = {SA)t + 
{SB)t- The z-component relaxes as a sum of two exponentially decaying quantities 
with different rates (corresponding to A and B). Therefore we cannot associate to it a 
total a single decay rate. 

The total transverse magnetization is the sum of that of species A and B. Each 
contribution evolves according to the modified Bloch equation. For large times, the 
dcphasing time approaches a renormalized constant value. Being again a sum of two 
terms decaying at different rates, the total transverse magnetization does not have a 
single decay rate. 

4 Proof of Theorem 2.1 

We set J = 1 in this proof (the case of general j is obtained merely by a change in 
notation). Following the method developed in [16], the dynamics of Ai is represented 
as 

(Ai), = (Vo, 5i • • • 5jve"^^if^5 n^^) . (4.1) 
The scalar product on the r.h.s. is that of the GNS Hilbcrt space ("doubled space"). 
Here, = ^Si <Si - • •<Si ' — ^Ri '^S' - • - f^Rjv ^^R and is the trace state of Sj , 
(fiSj ) ^^Sj) = 5([^]ii + [^]22), and Or^. are reservoir equilibrium states at temperature 
T = 1/13. 

The Bj are unique operators (in the commutant of the algebra of observables of 
spin j) satisfying 

^S,=Bjns„ (4.2) 

where is the initial state of Sj. 

The operator K is the Liouville operator acting on all spins and all reservoirs, 
satisfying 

Kng(g>n^ = o. (4.3) 

Its explicit form is easily written down (even though it is somewhat lengthy, see [16]) 
The main property is the representation 

Pr^'^Pr = + 0(a^e-^*), (4.4) 

e,s 

where = |Q^)(r2^| projects out all degrees of freedom of the reservoirs. The sum 
runs over all e of the form (2.8), i.e., eigenvalues of the operator 

Lg = Hg®lg- Ig® (4.5) 
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acting on C C (which are also the eigenvalues of K with a = 0). For each e 
fixed, s indexes its splitting into 1 < s < mult(e), as an eigenvalue of K, under 



the perturbation (2.2) plus (2.3).'^ We have ei^^ 7^ e'f ■* unless s = s'. The Qe^' are the 
(not orthogonal) spectral projections of K, and 7 > satisfies < ^ei^^ < 27 < T 
(temperature). 

(si (s) 

We now describe the perturbation expansion in a of eg and Ql . Due to As- 
sumption (2.9) the eigenspace of associated to an eigenvalue e is obtained as fol- 
lows. Associated to e are unique indices 1 < Ji < ^'2 < • • • < jNo(e) ^ ^ (recall 
(2.12)) satisfying aj = tj <^ j e {ji, ■ ■ ■ JNoie)} for any {a,f) with e{a,T) = e. Let 
Q = {qi, ■ ■ ■ 1 QNo{e)) € -|-l}-^°(^) and define vectors in by 



= 

= If 



a 



± 

3 J 



l + c,-[(af)*]2 



(4.6) 
(4.7) 



according to whether Qj = ±1. Here, Cj is given in (2.16) and 



with aj, bj, zf from (2.28), (2.15) and (2.26). 
Given e, ^, set 



(4.e 



'/e 



fai ,ri 



> <S) ■ ■ ■ ^ 



(?JVo(e) 
'i]Vo(e) 



where ^au,Th — Va-^. ® fr^ ^ C^, and where at locations jk, we replace ip, 
^) with the appropriate value of g^. 
Let hhe a form factor. We define 

Qh{u) = j |/i('U, S)pdS, and 7+(/i) = lim uQh{u). 

Jg2 u^Q+ 



(4.9) 
(4.10) 

by ^ (or 
(4.11) 



Let {/in} and {q!„} be form factors and coupling constants, respectively. For an eigen- 



*To be more precise, one has to use a 'spectral deformation' Kg of the operator K in this argument 
[16], but the deformation does not influence the physical results. 
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value e as in (2.8), set 

Xe({«n},{/in}) = T, «n^n P-V. / n^^^^ COth(/3H/2)du (4.12) 

{nrCTnT^Tn} 



ye({an},{/in}) = g XI «n('^n)^ ^/i„(wn) coth(/3a;„/2), (4.13) 

{n:(Tn^T„} 

^ ^ ^^^7+(/n), (4.14) 
{n:cT„^T„} 

^7+(/c)[eo(e)]^ (4.15) 



^o(e) = XI ^n(crn-Tn). (4.16) 
{n:(Tn^r„} 

Note that the indices over which the sums are taken are the same for any pair of spin 
configurations (£,r) with e(a, r) = e. Furthermore, we define 

Xe = Xe{{Xn},gc)+X,X{fin},{9n}), (4.17) 

Ye = ye+ye + ye{{K},9c)+ye{M,{9n}). (4.18) 

Then we have: 

Proposition 4.1 Suppose that the numbers e + si-, where 

No{e) 

6i^^=Xe + iYe+ X (4-19) 

k=l 

are distinct for all e and all g (the z are given in (2.26)). Then, for nonzero, small a, 
the eigenvalues of (the spectrally deformed) K are all simple and have the expansion 

e^i^ = e + <5?^ + 0{a^) (4.20) 

with corresponding eigenprojection 

g?^ = |r/?^)(#|+0(a^). (4.21) 



We give a proof of the proposition in Section 4.1. Combining the result of the Propo- 
sition with (4.1) and (4.4) gives 

(^i)t = E E (^Si • • • Vs., 5i • • • Bjv(|r/?^)(# 1)^1%) + Oia% 

e ^e{±l}^oM 

(4.22) 

with a remainder term uniformly bounded in t > 0. Since r^- belongs to the range of 
the spectral projection P{Lg = e), and since 

AiQg = P{Ls, = • • • = Ls. = 0)Ai%, (4.23) 
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only the terms e G spec(Li) = {— wi, 0, 0, CJi} in the sum in (4.22) contribute. 
Let us first consider e = — wi. We have A'^o('^i) = N — 1, 

r?il=<^+_®er<^---<^C and = <^+- ® ® • • • ® C'. (4.24) 

The term with e = —uji in (4.22) equals 

^it[_a,i+x_^i+iy_^,+Ef=2 4'+0(a^)] (4.25) 

N ^ 

X [p?^]2i \AM n (^s, , B,^ ) (ef , r2s, ) (4.26) 

.7=2 

= e^*(-i+^i+i^i+0(a^))c^(^^i) [pf ^]2i[Ai]i2, (4.27) 

where we set 

Ci(7V,t) = J] [e^*^; (^s„S,e+) (c;,51s,) +e"^^" (V's,,^,-^-) (5"'^s,)] • (4.28) 

Let us analyze the factors of this product. Using (4.6) and (4.7) we have (omitting the 
index f) 

{^S,Be)(e,^s) (4.29) 

= ^ _^ l[ae]'^ (V'S, -8(^11 + a^</?22)) (^(pii + c[q;^] V22, 2~^/'^{ipii + (^922) 

1 + ca^ I „ ^ 1 /o , 
1 + c[a^J^ \ 

(V;s,S{|(^i)((^i|®]l + a^|yj2)(<^2|<8l}^^s) (4.30) 

([po]ii + a^[po]22). (4.31) 



1 + c[Qt-']2 

1 + ca^ 
1 + c[ae]2 



In the last step, we use that B commutes with <8) 1, and that = V'S- Next 

we note the relation Q^a~ = — 1/c, which can be derived readily, for instance from the 
fact that + |^~)(^~| = 1. A short calculation then shows that 

^ 1 + ca+ , 1 + ca~ 



l + c(a+)2 l + c(a-)2' 

so that the factor in the product (4.28) becomes 

[po]ii {e"^\ + e'*^" (1 - 0} + (1 - [po]ii) {e'*^^a+C + e^*^"a-(l - C)} • 
Next, collecting the terms proportional to [po]ii and using 

(l-a-)(l-C) = -(l-a+)C and a+C = 1 - ^-(1 - C), 
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we obtain formula (2.25). 

One can transfer the error term down from the exponent: with D = 0{a^) we have 

^-tY,+w _ g-tn ^ g-in ^^^^ {wr ^^^^^ 

|g-m+ti? _ I < e-tri ^^t\D\ _ -^j < e-*^ii|L>|e*l^l (4.32) 

(mean value theorem). Now for \D\ < Ca^ and Yi > co? > 0, the r.h.s. can be 
bounded from above as follows: let e > 0, then for a < ce/C, an upper bound is 

Cto4e-*"'^(i-^) < Ca2e-*"'^(i-2^)supxe-^^^ = ^e-*"'^^!-^^). 

x>o ece 

This gives that if Yi > cc? > 0, then for all e > and a small enough, 

The remainder depends on e. Taking e = 1/2 and a < c/(2C), we get 

e-*[^^+«("')l =e-*^i+0(a2) (4.34) 

(uniformly in t > 0). This gives the contribution (2.22). 

Similarly to (4.27), one shows that the term in (4.22) with e = coi equals (2.23). To 
derive this, one checks that under the change e i— )■ — e, the exponent in (4.27) undergoes 
a complex conjugation and a sign change, and Ci turns into its complex conjugate. 

Next we conisder e = in (4.22). We have A^o(O) = N and obtain two contributions: 
one associated with eq = 0, rjo = (see (4.3)), ?7o = ^i" (8) • • • ® and another 
contribution with £q ^'^'^'"'^ and 

^(-1,1,1,...) ^^_^^+^___^^+ ~<-i,i,i,...) ^^_^^+^___^^+_ ^435^ 

It is easily seen that these two contributions are (2.20) and (2.21). This completes the 
proof of Theorem 2.1 ■ 

4.1 Level shift operators and proof of Proposition 4.1 

The total Liouville operator has the form 

K = Lo + W, (4.36) 

where Lq = + is the free (non interacting) Liouville operator, and W contains 
all interactions and is of O(a^) (for more detail, see Appendix B of [16], Resonance 
theory of decoherence and thermalization). To every eigenvalue e of Lq we associate 
the level shift operator 

Ae = -PeWPe (Lq - e + iO+y^PeWPe, (4.37) 

where Pg is the spectral projection of Lq onto the eigenvalue e. The eigenvalues of Ag 
are the second order (in a) corrections to the eigenvalues of K under the analytic (in 



19 



a) perturbation W oi Lq, see also Section 5 of [16], Resonance theory of decoherence 

and thermalization. Moreover, if ei~'' = e + S(r^ +0{a^) are the eigenvahies bifurcating 
out of e for a ^ then the corresponding eigenprojections of K, to lowest order in a, 

are given by \r]i-){r]^-\, where 



is) 



M is) 



[Ae\ rje- = [Oe- \ r/e" 



(4.38) 
(4.39) 



and (rji- jffe''^ = 1- For more information on these facts, we refer to Section 6 of [16], 
Resonance theory of decoherence and thermalization. We are assuming here that all 

energies si- are different, so that the corresponding eigenspaces are one-dimensional. 
This is generically true in applications, but it is not necessary for our strategy to work, 
see e.g. Appendix A of [16], Dynamics of collective decoherence and thermalization. 

Below, we give the explicit form of the level shift operators associated with all 
eigenvalues e (given by (2. 8)). Each level shift operators splits into a sum 



A, = Af^ + Ag 



loc 



(4.40) 



of two operators associated with the local and the collective interactions. We find the 

spectrum and eigenvalues of the level shift operators. In view of the explanations given 
at the beginning of this section, this gives a proof of Proposition 4.1. 



4.1.1 Collective level shift operator 

Let e be an eigenvalue (2.8) and define 



= 7>fneo(e)P.V. / 
4 Jk 



\m\' 

R3 \k\ 



d^A; 



(4.41) 



Proposition 4.2 (Collective LSO) The collective level shift operator associated to 
e is given by 



Af" = Xe{{Xn}, 9c) + i[y': + ye{{Xn}, 9c)] + ^ ^colb 

{n:a„=Tn} 

where M^jj acts on span{(^_|_-|_, } (doubled Hilbert space of n-th spin) as 



(4.42) 



1 -1 




' 1 ' 




- rn 


0-1 



(4.43) 



A proof of this proposition is obtained along the lines of Proposition 3.7 of [16] (Dy- 
namics of collective decoherence and thermalization). 
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4.1.2 Local level shift operator 

Proposition 4.3 (Local LSO) The local level shift operator associated to e is given 
by 

Ai- = Xe({Mn}, {gn}) + Av'e + ye({Mn}, K})] + ^1-' (4-44) 

{n:(Tn=Tn} 

where M^^^ acts on span{(^+_|_, } (doubled Hilbert space of n-th spin) as 



1 



-1 



(4.45) 



A proof of this proposition is obtained along the lines of Proposition 5.1 of [16] (Reso- 
nance theory of decoherence and thermalization). 

Remark. The contributions to local and collective level shift operators coming 
from the energy-exchange interactions are the same. 



4.1.3 Proof of Proposition 4.1 



The explicit forms of Ag° and Ag given in Propositions 4.2 and 4.3, and relation 
(4.40) yield 



Ae = + i^e + ^ 
{n:CT„=T„} 



-1 1 



1 
-1 



(4.46) 



The results (4.20), (4.19) and (4.9) for the eigenvalues and eigenvectors follow, and 
similarly for its adjoint. 

Note that for e = we have eo(e) = 0, and the level shift operator becomes 



TV 



Ao = i ^ 6„ 



n=l 



-1 1 



It follows that = and z„ = i6n(c„ + 1), with 



1 



V2 L 1 



-1 



V2 



1 



1 
-1 



(4.47) 



(4.48) 



5 Validity of perturbation expansion 

The Heisenberg equations of motion corresponding to (2.1)-(2.3) are 

5^ = -i[S'^,H] = X,Sy®ct>,{g,)+tinSy®M9n) (5.1) 

= UnSy-XnSy®Mfc)-'^nS^®Mfn) (5.2) 

Sy = -U;nS^-XnS^®M9c) + ^nS^®Mfc)-l^nS'^®M9n) (5-3) 

+l^nSn^Mfn)- (5.4) 
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For the local and collective annihilation operators a„(A;) and ac{k) we have 

ian{k) = [anik),H] = \k\an(,k) + ^iingn{k)S^ + ^Unfn{k)S^ (5.5) 

N N 

id,ik) = \k\a,{k) + -^Xngc{k)J2Sn + -^^nUk)Y,S'n- (5-6) 
^ n=l ^ n=l 

The latter two equations can be integrated, 

anik,t) = e-'^''Kn{k,0) 

-i= £ e-i|^l(*--) {Mn5n(A^)5^(T) + ,ynfn{k)S^nir)} dr (5.7) 
ac(A:,t) = e-'l*^l*ac(A;,0) 

"Tfi) ^"''"^*"'M^-5c(fc)i;5^(r) + x„/,(fc)^5^(r)|dr. (5.8) 

Remembering that 0(/i) = J^3{h{k)a* {k) + h{k)a{k)}d^k we insert (5.7), (5.8) into 
(5.2) to obtain 

+\>CnSl{t) [ /^-l^l(*--)^(A;)|A„5c(A;)f;S^(r) + x„/,(A;)f;S^(r)ldrd3A; 



2 

+h.c. 



+l'^nSm I f e-'l'l(*-"Vn(A^) {/Xn5n(A:)5^(r) + uMk)S'n{r)} drd^/c 

+h.c. (5.9) 

Let us denote by h'^ik) any of the products of functions of k occuring in the above 
integrals (e.g. h?'{k) = fc{k)gc{k) etc). Let us analyze the /c-integrals in the last 
expression for S"^. The product of form factors behaves like /i^(fc) = |A;|Pi+P2g-|'=l/*;o^ 
where pj = —1/2 + nj, rij = 0, 1, 2, . . . (see Assumption (C)). So 



[ ^^''^^h^{k)d^k ~ r rl+^i+'^^eir-rg-r-Aod^ = ^l+ni+n2 M 

Jr3 Jo 1/5 



l+ni+n2 



ro + ir 

which decays at least as (worst case ni = n2 = 0). Together with the boundedness 
|5'n'^'^(T)| < 1/2 this implies that the integrals over r and k in (5.9) are bounded 
homogeneously in i > 0, leading to 

For the validity of perturbation theory homogeneously in t > 0, we should impose 
x^iV, >CnXnN « ujn and jy^, VnfJ-n ^ ^n- Denoting by Uc and ac the size of collective 
and local coupling parameters, we thus need a^N « w, « u, where uo is the 
(typical) Bohr frequency of the single spin. 
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